Influence of the disorder on tracer dispersion in a flow channel 
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Tracer dispersion is studied experimentally in periodic or disordered arrays of beads in a capillary 
tube. Dispersion is measured from light absorption variations near the outlet following a steplike 
injection of dye at the inlet. Visualizations using dye and pure glycerol are also performed in similar 
geometries. Taylor dispersion is dominant both in an empty tube and for a periodic array of beads: 
the dispersivity Id increases with the Peclet number Pe respectively as Pe and Pe ' 82 and is larger 
by a factor of 8 in the second case. In a disordered packing of smaller beads (1/3 of the tube 
diameter) geometrical dispersion associated to the disorder of the flow field is dominant with a 
constant value of Id reached at high Peclet numbers. The minimum dispersivity is slightly higher 
than in homogeneous nonconsolidated packings of small grains, likely due heterogeneities resulting 
from wall effects. In a disordered packing with the same beads as in the periodic configuration, Id 
is up to 20 times lower than in the latter and varies as Pe a with a = 0.5 or = 0.69 (depending on 
the fluid viscosity). A simple model accounting for this latter result is suggested. 

PACS numbers: 47.56.+r,05.60.Cd 

I. INTRODUCTION are considered and a rich variety of dispersion behaviors 

is observed. 

A. Objectives of the paper 



Understanding mass and solute transfer in porous and 
fractured media is relevant to many fields of science and 
engineering due to their applications in domains such as 
waste managament and hydrology 0, 0, 0, • A partic- 
ularly sensitive method for characterizing mass transfer 
and detecting flow heterogeneities is the dispersion of a 
tracer in a flow of fluid through these media. 

In homogeneous 31? porous media, the variation in 
the flow direction x of the concentration C of passive 
tracer is often observed to satisfy at a macroscopic scale 
(a few pore sizes) the classical convection-diffusion equa- 
tion fl: 
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Here, U is the flow velocity and JDy is the dispersion co- 
efficient measured in the flow direction and which char- 
acterizes the longitudinal spreading of the tracers during 
their transport. 

In porous media, the dispersion coefficient is observed 
to be velocity dependent Q. At a low velocity, molec- 
ular diffusion predominates and the dispersion coeffi- 
cient is close to the molecular diffusion coefficient D m . 
When the velocity increases, different regimes are ob- 
served and non trivial relations such as power law vari- 
ations of the dispersion coefficient with the velocity are 
reported |1 S H H E3 • The objective of the present 
paper is to report dispersion experiments on well con- 
trolled models consisting of beads placed inside a long 
capillary. Different ordered and disordered bead layouts 



B. Key dispersion mechanisms 

Two physical mechanisms cause the dispersion of pas- 
sive tracers: molecular diffusion across and along stream- 
lines and advection. Advection refers to tracer transport 
by the motion of the host fluid. In order to compare 
both mechanisms, it is useful to define the Peclet num- 
ber: Pe — Ud/D m where d is a characteristic length 
scale. In this work, d is either the bead or the tube di- 
ameter (the latter is used when the flow tube does not 
contain any beads). 

In porous media, the void structure is complex, re- 
sulting in tortuous streamlines along which the velocity 
may greatly fluctuate. In randomly packed beads, for in- 
stance, current tubes often get split so that the fluid may 
flow around a bead; this results in variations of both the 
magnitude and the direction of the local velocity [TlllT^ |. 
With respect to the average flow, tracers therefore per- 
form a random walk between pores characterized by the 
duration t ~ d/U of single steps where d is the bead 
diameter. Using this description, one expects a longitu- 
dinal dispersion coefficient Dn/D m ~ d 2 /D m t ~ Pe and 
a dispersivity, l& = Du/U, constant and corresponding to 
the correlation length of the velocity field. Yet, regions of 
slow moving fluid, such as boundary layers, may also play 
a significant part and lead to logarithmic corrections |l3j| . 

In some specific configurations such as flow in a 
capillary tube (13. ITHl , in a periodic square array of 
beads or in a fracture [HIH0) the velocity field has 

long range correlations. In such media, the flow velocity 
difference between the walls bounding the open space and 
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its center stretches the tracer front and creates concen- 
tration gradients: The latter are balanced by molecular 
diffusion across the gradient. In such a case, the typ- 
ical time for the decorrelation of the tracer velocity is 
: t ~ d 2 /D m i.e. the characteristic time of molecular 
diffusion across the streamlines and the dispersion coef- 
ficient scales like Dn/D m oc Pe 2 . This regime is often 
called Taylor dispersion and, for flow in a capillary tube 
of inner diameter Z?, one finds |15| : 
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where Pe = UD/D rn . Equivalently, one can write : 
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While periodic 31? porous media are very difficult to re- 
alize experimentally (even using perfectly monodisperse 
beads), the transition from Taylor to geometrical disper- 
sion could be demonstrated numerically on 2D networks 
of increasing disorder |10| . 

Yet, some experiments on geological materials, either 
at the lab or field scales, display non Fickian charac- 
teristics. In these situations, the breakthrough curves 
observed are characterized by tails at long times. This 
effect may be explain either by the presence of hetero- 
geneities at the scale of the sample or by the trapping of 
solute in microscopic dead ends or stagnant pores. 

In this work, dispersion experiments in bead pack- 
ings of varying degrees of disorder are reported. De- 
pending on the structure of the packings, a variety of 
dispersion regimes is observed ranging from geometrical 
to Taylor dispersion and including intermediate regimes 
and non Fickian dispersion. From these experiments, the 
tube/bead diameter ratio appears as a key control pa- 
rameter of the dispersion regime. 



is D = 3.1 ± 0.1 mm, its length is 1500 mm and well 
calibrated stainless steel spheres of diameter d are used. 
Beads of diameter d = 2.54 ± 0.02 mm were chosen in 
order to reproduce the configuration used by Baudet et 
al [I?) . In addition to these two bead packs, two more 
configurations were used (see Fig. . The first one uses 
the same tube but empty. In the last one, the diameter of 
the beads is d = 1±0.01 mm, i.e. three time smaller than 
the tube diameter, leading to a larger tube-to-particle ra- 
tio. In order to clarify the present paper, the following 
convention is used : 

• (E) Empty channel: the capillary tube is devoid of 
beads. 

• (O) Ordered channel: the tube contains a periodic 
line of beads touching each other. 

• (DI) Disordered channel P. the beads are still 
touching each other but build up a disordered ar- 
ray. The tube/bead diameter ratio is 1.22. 

• {DII) Disordered channel IP. the arrangement is 
of the same type as (DI) but the diameter of the 
beads is smaller d = 1 ± 0.01 mm. The tube/bead 
diameter ratio is 3.1. 

In configurations (O) and (DI), the diameter of the 
spheres is slightly smaller than the tube diameter (d — 
2.54 ± 0.02 mm). In configuration (DII), the diameter 
of the beads is three times smaller than that of the tube 
(d = 1 ± 0.01 mm) Originally, one end of the capillary 



PROCEDURE 




op 






VP 





















II. 



The experimental model used in this work closely re- 
sembles that used by Baudet and coworkers 01- ^ n this 
latter work, tracer dispersion was measured in a long 
capillary tube filled with monodisperse beads and with a 
tube/particle diameter ratio of the order of 1.25. The lay- 
out of the beads was either periodic or disordered. The 
first configuration was achieved by adding one by one 
beads in an horizontal capillary tube. Finally, one ends 
up with a line of beads, each of them being in contact 
with its two neighbors and with the tube wall. The sec- 
ond layout is obtained by tilting the capillary tubes while 
filling them. The beads fall then on top of each other and 
slide often sideways, resulting in a denser packing. The 
final packing is random and beads are added in the tube 
until it is entirely filled. 

In this work, the inner diameter of the capillary tube 



DI 




FIG. 1: Flow geometries used in the experiments. From top 
to bottom : Empty Channel - Ordered Channel - Disordered 
Channel I - Disordered Channel II. 

tube is connected to a syringe pump allowing to estab- 
lish a stationary flow of the transparent fluid through the 
channel: A valve allows to switch the injection to a dyed 
solution of identical properties. Fluid flowing out of the 
other end of the tube is weighed by computer controlled 
scales, allowing to measure the flow rate throughout the 
experiment. 
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The fluids used are solutions of either 10% or 70% of 
glycerol in water and their dynamic viscosities /i are re- 
spectively 1.3 x 10~ 3 Pa.s and 23 x 10~ 3 Pa.s (at 20°C). 
Water Blue dye at a concentration of 0.05 g/1 is added 
to one of the solutions: It has been selected because it 
is chemically stable and does not modify the rheologi- 
cal properties of the solution. The molecular diffusion 
coefficient D m of the dye was determined through inde- 
pendent Taylor dispersion measurements performed in a 
vertical teflon capillary tube. One obtains in this way 
D m = 6.5 10~ 4 mm 2 /s for the 10% glycerol solution; 
for the 70% solution, the value D m — 3.25 10~ 5 mm 2 /s 
is deduced from the value for the first solution by the 
Stokes- Einstein relation [l3 |. In the present experiments, 
the Peclet number ranges from 50 to 10 4 for the 10% 
glycerol solution but reaches 1.4 x 10 5 for 70% glycerol 
solutions. In all the studies the Reynolds number is less 
than 1 except for the empty tube where Re can reach 
10: Under such conditions, the flow can be considered as 
stationary |ll| . 

In this work, the variation of the tracer concentra- 
tion at the outlet with time, known as the breakthrough 
curve, is determined from light absorption by the dye. 
The measurement is realized over a square window of size 
0.3 mm 2 located at 1450 mm from the injection, close to 
the outlet and immediately downstream of the last bead 
of the bed. In order to reduce optical distorsion induced 
by the tube curvature, the measurement section is en- 
closed within a transparent plexiglas cell with flat paral- 
lel walls. The cell is originally filled with glycerol, a fluid 
with a refractive index close to that of the tube. The 
section of the tube inside this cell is inserted between a 
light panel and a 4096 gray levels CCD camera (Roper 
Coolsnap CF). The set-up is illuminated by a fluorescent 
tube placed on the opposite side from the camera and 
excited at a high frequency to reduce fluctuations. For 
each experiment, 2000 images are recorded by a com- 
puter connected to the camera at time intervals ranging 
from 1 to 30 sec depending on the flow rate. Dye concen- 
tration values are determined quantitatively using cal- 
ibration measurements realized independently with the 
experimental tube saturated with dye solutions of dif- 
ferent known concentrations. Finally, drifts of the light 
intensity are measured in a region of interest outside the 
tube; these measurements are then used during the anal- 
ysis of the images to compensate for the effect of these 
variations on the transmitted light intensity in the exper- 
imental section. 

Using this experimental procedure, breakthrough 
curves were measured for different flow velocities and for 
different bead layouts. The results of these global mea- 
surements are given and discussed in section llTTI More- 
over, in order to improve interpretations of the break- 
through curves, visualizations realized independently at 
a local scale in a similar configuration were carried out. 
In these latter experiments, a transparent tube of D = 
8 mm inner diameter is filled with glass beads of diame- 
ters 6, 3 or 2 mm. These beads were chosen so that the 



corresponding tube/bead diameter ratio is respectively 
1.33, 2.67 and 4, close to the values used in the disper- 
sion experiment. Originally, the model is saturated with 
glycerol (viscosity ~ 1 Pa.s) which is displaced by the 
same fluid but dyed. Using a fluid with such a high vis- 
cosity (and therefore a low diffusion coefficient) allows to 
observe and separate clearly the various flow paths at low 
Reynolds numbers inside the sample. 

Before describing in detail the various dispersive 
regimes observed for the different bead layouts, the next 
section describes the methods used to determine the dis- 
persion coefficients from the breakthrough curves. 



A. Analysis of experimental curves 

Figure displays a typical breakthrough curve ob- 
tained after a stepwise injection of the dyed fluid; asym- 
metrical curves were also obtained and their analyzis 
will be discussed later. Under such initial conditions 
and if the concentration satisfies the classical convection- 
diffusion equation given by Eq. the concentration 
variation at the outlet is given by: 



c 2 U J ^m^t h 



(4) 



Where, Co is the dye concentration in the displacing fluid, 
L the distance between the measurement section and the 
injection point, Dm is the longitudinal dispersion coeffi- 
cient, U is the mean flow velocity. The figure J2Jl shows 
the fit of the breakthrough curve by the function given 
by Eq. |0J where the only adjustable parameter is D\\. 
The two curves are almost undistinguishable, indicating 
that the dispersion process is Fickian and that the clas- 
sical convection-diffusion equation is satisfied. Yet, 
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FIG. 2: Continuous line : typical tracer variation concen- 
tration as a function of time for a type DI— array and for 
Pe — 930 (continuous line).- Dotted line : curve fitted with a 
variation from Eq. Jlj) 

such perfect fits were not always observed. For instance, 
in the O configuration and for fairly high velocities U, 
one observes at long times, as can be seen on Fig. [21 a 
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" tail" effect : it corresponds to an exponential relaxation 
of the concentration towards a limiting value with a typ- 
ical time of the order of L/U. Such breakthrough curves 
are better fitted by solutions of the classical Coats-Smith 
capacitive model |l9j which uses four fitting parameters: 
t = L/U, D«, f and tf. The two last parameters cor- 
respond to the amplitude and the characteristic time of 
the exponential variation in the "tail" at long times. By 
combining these parameters, one obtains an asymptotic 
dispersion coefficient -D|| as : 
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(5) 



D\\ as represents the value of the dispersion coefficient 
that would be measured for channels with the same local 
structure, but long enough so that a Gaussian dispersion 
regime is reached and Eq. becomes valid. Note that 
the Coats-Smith model is only used here as a mathemat- 
ical way of obtaining the value of D\\ as ; the fact that the 
experimental curves are well fitted by the model does not 
mean that its underlying assumptions are valid (namely 
the existence of a zero flow zone with an exponential ex- 
change with the mean flow). In the following, the disper- 
sivity Id is taken equal to Du as /U only for experiments 
in the O configuration, in the other configurations, D\\ is 
obtained by fitting the experimental curves to Eq. 
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FIG. 3: Continuous line : experimental dye concentration 
variation as a function of time for an ordered O— type array 
of 2.54 mm diameter beads and for Pe = 1000. Dotted lines 
: Gaussian and Coats-Smith fits. 



III. EXPERIMENTAL RESULTS 

Figure 0] displays variations of the dispersivity Id as a 
function of the Peclet number measured for the different 
channel configurations used in the present work. Note 
that all the experiments were carried out twice: First, 
the dyed fluid was injected to displace the clear solution 
flowing initially in the model at the same flow rate. After 
the model has been completeley saturated with the dyed 
fluid, it is displaced in a second experiment by the clear 
fluid, still at the same flow rate. The dispersion coeffi- 



cient was found to be identical in the two experiments, 
implying that no instabilities modified the flow. 
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FIG. 4: Variation of the dispersivity Id as a function of the 
Peclet number: (A) DII— array; (o, •) DI— array respectively 
for solutions of 10% and 70% glycerol in water; (■) E— array; 
(□) O— array (in this latter case Id = D\\ a3 /U). 



A. Empty tube (E— Channel) 

In the case of an empty tube, all the experimental 
dispersion curves are well adjusted by the concentration 
variation given by Eq. The corresponding disper- 

sivities Id, as plotted in Figure 01 increase linearly with 
Pe reflecting, as expected from Eq. Q, a dominant in- 
fluence of Taylor dispersion. One sees from Eq. J3J| that 
the value of Id is related to the tube diameter; a linear 
regression of the data gives indeed an effective value of 
the tube diameter : a e ff — 3 ± 0.2 mm. 



B. Ordered array of beads (O Channel) 

Dispersivity Id values for the ordered array of beads 
plotted in Figure 0] were obtained with the Coats-Smith 
model. For Pe < 2500, Id varies with Pe following a 
power law : Id ~ a e ffPe a (straight line in log-log coor- 
dinates). From a regression over the experimental data, 
one obtains a e ff — 0.6 mm and a = 0.85 ± 0.2: This 
value is compatible with a Taylor dispersion mechanism 
for which a = 1. 

One can indeed expect that, in this channel, fluid 
particles follow streamlines determined by the periodic 
structure of the bed of beads: Therefore, as in capil- 
lary tubes, tracer can only move from a streamline to 
another through molecular diffusion (See Fig. [5J . An im- 
portant feature is the fact that the experimental values 
of Id which vary from 40 mm to 400 mm are about 8 
times larger than for the empty capillary tube at a same 
mean flow velocity. This cannot be explained by the par- 
tial filling up of the flow channel by the beads. On the 
opposite, since these reduce the effective aperture of the 
flow channels, Id should drop off following Eq.© which 
is not observed. 
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Direct visualizations realized in a similar periodic ge- 
ometry (Figure [SJ, although with larger beads, help un- 
derstand these results: Pure glycerol is used in these ex- 
periments to visualize clearly the boundaries between flu- 
ids by removing the influence of molecular diffusion. In 




FIG. 5: Miscible displacement of pure by dyed glycerol in an 
ordered array of beads inside a capillary tube : (a) side view 
(b) top view (bead diameter : d = 6 mm; tube diameter : 
D = 8 mm) 

this figure, one observes that the dyed fluid displays a 
tongue like structure and wraps around the beads. This 
shape shows that the streamlines are mostly oriented in 
the flow direction with tiny undulations induced by the 
beads. Another important point is that the dyed fluid 
does not flow in the narrow space between the beads: 
Mass transfert between the two regions can thus only oc- 
cur through molecular diffusion. This explains the long 
tail observed in the breakthrough curves. Finally, the 
tongue-like structure reveals the strong velocity contrast 
between the fluid flow on the side of the beads and above 
them. The velocity gradient stretches the tracer fronts 
resulting in a concentration gradient which, in turn, gets 
smoothed by a transverse diffusive flux. 
All ingredients of a Taylor like dispersion regime are thus 
present. Yet, because of the beads, the diffusive flux is 
more tortuous than in an empty capillary tube: Tracer 
has to flow around the beads to reach the outlet of the 
tube. The time needed to homogenize the tracer concen- 
tration in the tube section is thus longer, resulting in a 
higher dispersivity. 

At the highest flow velocities, the variation of Id lev- 
els off and it starts to decrease. A possible explanation 
is the fact that, in this range of Pe values, the Reynolds 
number becomes higher than 1 (Re > 5): recirculation 
flows may then develop and induce a more efficient trans- 
verse mixing than molecular diffusion. This may account 
for the decreasing value of Id at high Peclet numbers. 



C. Disordered array of small beads (DII —array) 

For a disordered array of 1 mm diameter beads inside 
the tube, the dispersivity Id increases slowly as a function 
of Pe before reaching a constant value of the order of 
3 mm for Pe > 600. This constant limiting value of Id 
implies that geometrical dispersion associated with the 
random velocity variations from one pore to the next 
is dominant. Also, in such geometries, the correlation 



length of the velocities of fluid particles along their path 
is too short for Taylor dispersion to have a significant 
influence. Dispersion characteristics of such arrays are 
then comparable to those of non consolidated packings 
of monodisperse grains: In this latter case, the minimum 
value of Id is lower than the grain size (typ. 0.6 — 0.7 d g ) 
and is reached for Peclet numbers (based on the grain 
size) of the order of 10 20]. In Fig. QJ the value of Id 
for a similar Peclet number (Pe = 30 when based on the 
tube diameter) is Id = 2 mm or about twice the bead 
diameter. This result also suggests that boundary layers 
have a weaker influence on dispersion than predicted by 
some theories since no increase of Id as LogPe is 
observed. 




FIG. 6: Miscible displacements of pure glycerol by dyed 
glycerol saturating disordered glass beads packings inside a 
D = 8 mm diameter capillary tube. Bead diameter : (a) 
d — 2 mm - (b) d = 3 mm. The tube/particle diameter ratio 
is respectively 4 and 2.66. 



Visualizations realized with d = 2 and 3 mm beads 
inside a D — 8 mm diameter tube (FigsHJ^-b) help under- 
stand these results. Figure demonstrates clearly the 
many divisions of the front of injected fluid after mov- 
ing through several pores and the rather uniform dis- 
tribution of the invading fluid across the flow section : 
for the less viscous fluids used in our dispersion experi- 
ments, transverse molecular diffusion would mix quickly 
these thin filaments with the surrounding fluid, leading 
to the Gaussian dispersion observed for the DII— array. 

For the lower value 2.66 of the ratio d/D the shape of 
the dye streaks changes very much as shown on FigureHJjD. 
In this particular case, flow is strongly channelized in the 
center of the column; this indicates, that for some (low) 
values of the ratio d/D, pathological structures of the 
flow field may appear. The dispersion measurements re- 
ported in the present section were performed with a col- 
umn for which d/D was intermediate between the two 
values in Figs. HJ^-b. Since no channelization of the flow 
was observed visually during the measurement, the geo- 
metrical dispersion observed corresponds most likely to 
the type of flow structure shown in FigHJi. In the next 
section, tracer dispersion is investigated for still smaller 
values of d/D. 
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D. Disordered array of large beads (DI— array) 

In this case, the beads have the same diameter as in 
the periodic O— channel (i.e. d = 2.54 mm), but the 
packing is now disordered (See Fig.^l. The experiments 
were performed using the two water-glycerol solutions 
containing either 10% or 70% of glycerol in weight. As for 
the DII array and in contrast with the ordered one, all 
breakthrough curves are well adjusted by Eq. allowing 
for the determination of Du for the various flow condi- 
tions. The Figure 0] shows the dispersivity Id = D\\/U 
as a function of the Peclet number for the two solutions 
used in the experiments. 

Clearly, the dispersivity Id and Pe have a power law rela- 
tionship. Fitting the variation of Id with Pe to a e ffPe a 
gives respectively for the 10% and 70% glycerol solutions 
a = 0.52 ± 0.01, a eff ~ 0.2 mm and a = 0.69 ± 0.04, 
a e ff — 0.04 mm. One observes that a e // and a are only 
slightly different for the two solutions despite a ratio of 
20 between the viscosities. One notices that the values 
of Id for a given Peclet number are much lower for the 
disordered array than for the periodic one (by a factor 
of 15 (resp. 40 at low (resp. high) Pe values). Also, at 
low Peclet numbers, Id becomes of the order of the bead 
size and gets close to the dispersivity observed for the 
DII array. 

An important feature is the fact that the exponent a 
characterizing the variation of Id with Pe is of the order 
of 0.5 (for the 10% solution) and of 0.69 for the other one 
(at slightly higher Peclet numbers). These exponents are 
intermediate between the values 0.82±0.2 for the periodic 
array and 1 for the empty tube and the value (constant 
Id) for the disordered packing of smaller beads: This im- 
plies that the dispersion mechanism is intermediate be- 
tween Taylor and geometrical dispersion corresponding 
respectively to the first and second cases. This result 
was also observed by Baudet et al 01 but remained un- 
explained. 

As in the previous sections, visualization on a larger sys- 
tem with the same d/D ratio help understand the disper- 
sion mechanisms. The visualizations of FigureE^-b show 
that the dyed fluid is mostly split into two streaks located 
near the walls. This channelization results from the in- 
crease of the porosity near the walls |2^; moreover, the 
low value of the tube/bead diameter ratio (of the order 
of 1) breaks the angular isotropy of the porous struc- 
ture and concentrates the flow paths in a few (here 2) 
channels. Magnico ^2 estimated that, at low Reynolds 
numbers, a fluid layer of thickness of the order of d/4 ap- 
pears, inside which fluid flow is purely longitudinal and 
tangential with no radial component. At first, one ex- 
pects therefore radial exchange between the dye streaks, 
clearly visible on Figure^, and the remaining pore space 
to be mostly diffusional. Yet, as can be seen of the Fig- 
ure 03 structural heterogenety resulting from misplaced 
beads split the streaks into filaments parallel to the mean 
flow pl| ; the number of filaments increases then along 
the flow path and their size decreases. These filaments 



persist over distances significantly larger than the bead 
diameter so that molecular diffusion may spread the fila- 
ments over transverse distances of the order of their size. 
This gives rise to Taylor-like dispersion so that the global 
dispersion results from the combined influences of geo- 
metrical and Taylor dispersions. Such an influence of the 
flow channelization on dispersion was recently reported 
by Bruderer et al |10| in 2D networks. 




FIG. 7: Views at different times of the miscible displacement 
front of pure by dyed glycerol in a disordered channel: the 
diameter of the beads is d = 6 mm and the diameter of the 
tube D — 8 mm. The tube-to-particle ratio is 1.33. 



E. Qualitative model of different power law 
variations 

A qualitative argument helps understand how one can 
reach such power law variations of Id oc Pe a (0 < a < 1) 
under the combined influence of the disorder of the flow 
field and of transverse molecular diffusion. Assume that 
the front gets divided into streaks of width a x decreasing 
with the distance x parallel to the flow as a x ~ d 1+l3 3 



(the d term allows to have the right dimensionality for 
the equation). By generalizing the Taylor argument, the 
transition to diffusive spreading should occur when the 
transverse molecular diffusion time Tdiff across the dis- 
tance a x is of the order of the mean transit time L/U 
with : 



T diff 



D r . 



x 
U 



(6) 



Replacing a x by its expression provides the distance at 
which the transition should take place : 



x / ,. 



(7) 



With, as usual, Pe = Ud/D m . As in Taylor disper- 
sion, Xtransi represents the characteristic decorrelation 
distance of the velocity of tracer particles and we assume 
therefore that D» ~ Ux trans leading to 
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(9) 



7 



For (3 = (no geometrical variation of the filament 
width with distance), one retrieves Taylor dispersion with 
D\\Tayior ~ a 2 U 2 / D m and for (3 = oo (fast decorrela- 
tion), one obtains geometrical dispersion with D\\ geom ~ 
dU (the exponent a defined previously should then be 
related to (3 by a = 1/(1 + 2/3) In the present case, the 
experimental result oc Pe 3 / 2 implies that (3 — 1/2. 
Therefore, the experimental observations on DI— may 
be accounted for by assuming a combination of the influ- 
ences of transverse molecular diffusion and of the geomet- 
rical disorder of the packing with, for the latter, a rate of 
division of the fluid streaks intermediate between those 
observed in a packing of small beads and in a capillary 
tube. 



IV. CONCLUSION 

To conclude, despite its simple structure (a long tube 
filled up with beads), the experimental system studied 
in the present work displayed a broad variety dispersion 
regimes. Both the layout of the beads and the ratio 
of their diameter to that of the tube were shown to 
influence very much the dispersion characteristics. 
First, Taylor dispersion is observed for a periodic bead 
array inside a tube of slightly larger diameter than the 
beads, like in the empty tube but with an increased 
dispersion coefficient. With the same beads inside 
the same tube, but packed in a disordered array, the 
dispersivity Id is very strongly reduced and the exponent 
a characterizing the variation of Id with Pe decreases 



from almost 1 to 0.5. This variation reflects the reduced 
persistence length of fluid streaks as they move along 
the tube (this length remains however much larger than 
the bead diameter). A simple model assuming a power 
law reduction of the width of the streaks of dye with 
distance allows to reproduce this variation of Id with 
Pe. When the ratio of the tube and particle diameters 
is increased by using beads of smaller diameter, the 
dispersion coefficient varies with the Peclet number as 
D|| oc Pe for Pe > 600. This reflects a geometrical 
dispersion regime with a short persistence of the dye 
streaks controlled by the local pore geometry and a weak 
influence of diffusion in the boundary layers. 
More observations of the flow and concentration fields 
at the local scale (using for instance matched index 
fluids) are needed to explain quantitatively these results. 
The increase of the dispersivity for the periodic array 
compared to the empty tube and the decorrelation of the 
motions of the fluid particles in the disordered arrays for 
small tube/particle diameter ratios are two particularly 
important issues. 
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